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Abstract 

The azimuthal anisotropy of jet spectra due to energy loss of jet partons in azimuthally 
non-symmetric volume of dense quark-gluon matter is considered for semi-central nuclear 
interactions at collider energies. We develop the techniques for event-by-event analysing 
the jet azimuthal anisotropy using particle and energy elliptic flow, and suggest a method 
for calculation of coefficient of jet azimuthal anisotropy without reconstruction of nuclear 
reaction plane. 
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1 Introduction 



High-pr jet production and other hard processes are considered as the promising tools for 
studying properties of hot matter created in heavy ion collisions at RHIC and LHC The 
challenging problem is the behaviour of colour charge in quark-gluon environment associated 
with the coherence pattern of the medium-induced radiation. It results in interesting non-linear 
phenomena, in particular, the dependence of radiative energy loss per unit distance dE/dx 
along the total distance traversed (see review |]J and references therein). In a number of 
papers [fj, |3|, J|, |j authors analysed the azimuthal anisotropy of high-p T hadron spectra in semi- 
central nuclear collisions at RHIC due to partonic energy loss in azimuthally non-symmetric 
volume of quark-gluon plasma. At LHC energy, when the inclusive cross section for hard jet 
production on E T ~ 100 GeV scale is large enough to study the impact parameter dependence 
of such processes [[J, one can hope to observe the similar effect for hadronic jet itself J71. In 
particular, CMS experiment at LHC || will be able to provide jet reconstruction and adequate 
measurement of impact parameter of nuclear collision using calorimetric information || . 

It is important to notice that the coherent Landau-Pomeranchuk-Migdal radiation induces 
a strong dependence of the radiative energy loss of a jet (but not a leading parton) on the 
jet angular cone size [JTOj, [11], [12], [13|]. It means that the medium-induced radiation will, in 
the first place, soften particle energy distributions inside the jet, increase the multiplicity of 
secondary particles, and to a lesser degree affect the total jet energy. On the other hand, 
collisional energy loss turns out to be practically independent of jet cone size, because the bulk 
of "thermal" particles knocked out of the dense matter by elastic scatterings fly away in almost 
transverse direction relative to the jet axis [[TI[ . 



The methodical advantage of azimuthal jet observables is obvious: one needs to reconstruct 
only azimuthal position of jet, but not the total jet energy. It can be done more easily and with 
high accuracy, while the reconstruction of the jet energy is more ambiguous task ||. However 
the performance of inclusive analysis of jet production as a function of azimuthal angle requires 



event-by-event measurement of the reaction plane angle. The summarized in papers |LJ] present 
methods for determination of the reaction plane angle are applicable for studying anisotropic 
flow of soft and semi-hard particles in current heavy ion dedicated experiments at SPS [Tj| and 



RHIC JT6[, and, in principle, might be also used at LHC 0. 

In this work we suggest a method to calculate the coefficient of jet azimuthal anisotropy 
without reconstruction of nuclear reaction plane. In some sense, it represents the development 
and generalization of the well-known method for measuring azimuthal anisotropy of particle 
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flow considered originally in a number of works (see |n], [T7|, [TJj for instance). 



2 Event-by-event analysis of azimuthal correlations 

Let us remind the essence of techniques |17], [14]] for measuring azimuthal elliptic anisotropy of 



particle distribution, which can be written in the form 

% = |[l + 2, 2 cos2( V -«], K = f*p%. (1) 

Knowing the nuclear reaction plane angle i/># allows one to determine the coefficient v 2 °f 
azimuthal anisotropy of particle flow as an average (over particles) cosine of 2(p: 

1 } dN 

(cos2(ip - ip R )) = — J dip cos2(v9 -t/j R ) — = v 2 ■ (2) 

— 7T 

However in the case when there are no other correlations of particles except those due to flow 
(or such other correlations can be neglected), the coefficient of azimuthal anisotropy can be 
determined using two-particle azimuthal correlator without the event plane angle ip R : 

1 J J d 2 N 
(cos2(y3 1 - <p 2 )) = j-p J dip l J dip 2 cos2(y) 1 - cp 2 ) d d (3) 



dN dN 

= Vo 



— 7T — 7T 

Here one should to note that it is safe to neglect non-flow correlations only if the coefficient 
of azimuthal anisotropy v 2 is much larger than 1 / \/Nq. In reality, this condition is not always 
satisfied. A new method, based on a cumulant expansion of multi-particle azimuthal correla- 
tions, which allows measurements of much smaller values of azimuthal anisotropics, down to 



1/iVo, has been worked out recently in paper [19]. This method automatically eliminates the 



major systematic errors, which are due to azimuthal asymmetries in the detector acceptance, 
and in principle could be used in further improvements of our approach. 

Let us consider now the event with high-p^ jet (dijet) production, the distribution of jets 
over azimuthal angle relatively to the reaction plane being described well by the elliptic form [[?]], 

dN jet N jet j dN& 

= _ [1+2^ oob2( V -^)], Ni =jd<p — , (4) 

— 7T 

where the coefficient of jet azimuthal anisotropy v 2 et is determined as an average over all events 



cosine of 2ip, 



1 f , „, , x dN jet 



(cos2((p-ip R )) event = J dtp cos2((p-ip R ) ^— - = v\ et . (5) 

-77" 
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One can calculate the correlator between the azimuthal position of jet axis ipj e t f\ and the angles 
of particles, which are not incorporated in the jet(s). The value of this correlator is related to 
the elliptic coefficients v 2 and v J 2 et as 

1 } J dN jet dN 

«cos2(^ - ip))) event = -p- J df jet j dtp cas2(<p jet -ri-jj-fy (6) 

— 7T —TV 

1 } dN jet ■ 

dtp jet cos2(f jet -ifj R ) - v 2 = v° 2 v 2 ■ 



Using Eq. (^) and intermediate result in Eq. (Q) (after averaging over particles cos 2(fj et — ip) 
reduces to v 2 cos2(fj et — ipn)) we derive the formula for computing absolute value of the coef- 
ficient of jet azimuthal anisotropy (without reconstruction of sign of v 2 et ): 

v iet = I (cOS2(ip jet - ip)) \ ^ 

\ ^(cos 2{y 1 -y 2 ))l ewnt 

This formula does not require the direct determination of reaction plane angle ipR- The brackets 
( ) represent the averaging over particles (not incorporated in the jet) in a given event, while 
the brackets ( ) event the averaging over events. 

The formula (|7]) can be generalized by introducing as weights the particle momenta, 

y jet = I (cQS%jet ~<p) Pt(<P)) \ ^ 

\ v /(cos2(^i - <p 2 ) p T i(fi) PmM) I event 

In this case the brackets ( ) denote the averaging over angles and transverse momenta of 
particles. The other modification of @, 

jet I (cos2(^ jef -^) E(<p)) \ 
V 2(E) — \ r, ~ \ - ; — \ ~ ; — — / \ J ) 



(cos2(<^i - if 2 ) #1(^1) E 2 (f2)), 



event 



(Ei(fi) being energy deposition in calorimetric ring % of position ipt) allows one using calori- 
metric measurements (§) for the determination of jet azimuthal anisotropy, in particular, under 
condition of CMS experiment at LHC. 

3 Numerical simulation and discussion 

In order to illustrate the applicability of method presented for real physical situation, we con- 
sider the following model. 



lr rhe other possibility is to fix the azimuthal position of a leading particle in the jet. In this case calculating 
azimuthal correlations can provide the information about azimuthal anisotropy of high-pr particle spectrum. 



3 



Jets. The initial jet distributions in nucleon-nucleon sub-collisions at y/s = 5.5 TeV have 
been generated using PYTHA_5.7 generator pOfl . We simulated the rescattering and energy 
loss of jets in gluon-dominated plasma, created initially in nuclear overlap zone in Pb— Pb 
collisions at different impact parameters. For details of this model one can refer to our previous 
papers || fjj. Essentially for our consideration here is that in non-central collisions the azimuthal 
distribution of jets is approximated well by the elliptic form @. In the model the coefficient of 
jet azimuthal anisotropy increases almost linearly with the growth of impact parameter b and 
becomes maximum at b ~ 1.2Ra, where Ra is the nucleus radius. After that v J 2 et drops rapidly 
with increasing b: this is the domain of impact parameter values, where the effect of decreasing 
energy loss due to reducing effective transverse size of the dense zone and initial energy density 
of the medium is crucial and not compensated more by stronger non-symmetry of the volume). 
Other important feature is that the jet azimuthal anisotropy decreases with increasing jet 
energy, because the energy dependence of medium- induced loss is rather weak. Finally, the 
kinematical cuts on jet transverse energy and rapidity has been applied: E 3 ^ 1 > 100 GeV and 
\y ]et \ < 1.5. After this dijet event is superimposed on Pb— Pb event containing anisotropic flow. 



Particle flow. Anisotropic flow data at RHIC [JIB] can be described well by hydrodynamical 
models for semi-central collisions and pt up to ~ 2 GeV/c |2~lfl , while the majority of microscop- 
ical Monte-Carlo models underestimate flow effects (however, see ||22|| ). One can expect that 
the hydrodynamical model can be applied to estimation of particle flow effects at LHC, may 
be extending this approach to even higher p^ values. On the other hand, at collider energies 
one more reason for anisotropic flow in relatively high-p^ domain can arise: the sensitivity of 
semi-hard particles to the azimuthal asymmetry of reaction volume under the condition that 
the major part of them are the products of in-medium radiated gluons or parent hard par- 
tons P, [§. Of course, for more detailed simulation, one has to take into account the interplay 
between hydro flow and semi-hard particle flow. However here we restrict our consideration to 



using the simple hydrodynamical Monte-Carlo code |23| giving hadron (charged and neutral 
pion, kaon and proton) spectrum as a superposition of the thermal distribution and collective 
flow. For the fixed in the model "freeze-out" parameters — temperature 7/ = 140 MeV, col- 
lective longitudinal rapidity Y™ ax = 5 and collective transverse rapidity Y™ ax = 1 — we get 
average hadron transverse momentum < p\ >= 0.55 GeV/c. We set the Poisson multiplicity 
distribution and take into account the impact parameter dependence of multiplicity in a sim- 
ple way, just suggesting that the mean multiplicity of particles is proportional to the nuclear 
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overlap function (dN/dy) (b) oc Taa(P). In the framework of this model, anisotropic flow can 
be introduced on the assumption that the spatial ellipticity of "freeze-out" region, 

{y 2 + x 2 ) ' ^ 

is directly related to the initial spatial ellipticity of nuclear overlap zone, eo = 6/2i2^. Such 
"scaling" allows one to avoid using additional parameters and, at the same time, results in intro- 
ducing elliptic anisotropy of particle and energy flow due to dependence of effective transverse 
size of "freeze-out" region Rf(b) on azimuthal angle of "hadronic liquid" element <3>: 



Rf(b) = Rf(b = 0) min{ yl - e 2 , sin 2 $ + e cos$, yl - e 2 , sin 2 $ - e cos$}. (11) 

Obtained in such a way azimuthal distribution of particles is described well by the elliptic 
form (|J) for the domain of reasonable impact parameter values. Note that in the framework 
of present paper we do not aim at detailed description of azimuthal hadronic spectra and its 
comparison with experimental data (in principle, fit of the data could be performed using, for 
example, eo in ( |TTD as a parameter). In order to investigate the reliability of the method, we just 
need here to introduce the elliptic anisotropy of energy flow in Monte-Carlo event generator. 

Energy flow. To be specified, we consider the example of CMS detector at LHC collider ||. 
The central ("barrel") part of the CMS calorimetric system will cover the pseudo-rapidity 
region \r)\ < 1.5, the segmentation of electromagnetic and hadron calorimeters being At] x 
A<p = 0.0174 x 0.0174 and A77 x A<p = 0.0872 x 0.0872 respectively @. In order roughly to 
reproduce the real experimental situation (not including real detector effects, but just assuming 
calorimeter hermeticity), we apply formula @ to integrated over rapidity energy deposition 
Ei(ifi) of generated particles in 72 rings (according to the number of rings in the hadron 
calorimeter) covering full azimuth. 

Figure 1 shows the 6-dependence of " theoretical" value of v J 2 et (calculated including collisional 

and radiative energy loss), and v 3 2 et determined by the method (§) for the two estimated values of 

the input parameter, number of charged particles per unit rapidity at y = in central Pb— Pb 

collisions, dN±/dy = 3000 and 6000. One can see that the accuracy of v 2 e determination 

is close to 100% for semi-central collisions and becomes significantly worse in very peripheral 

2 In fact, the influence of initial nuclear effects like shadowing or saturation can result in a weaker impact 
parameter dependence of (dN/dy) (see for example |24|). Thus, if the particle density in central collisions is 
fixed as a parameter, we are getting here the most "pessimistic", minimal estimation of mean multiplicity of 
non-central events. Since the accuracy of the method used improves with the increasing multiplicity, for weaker 
^-dependences one can obtain even better resolution for v J 2 et . 
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collisions (b ~ 2Ra), when decreasing multiplicity and azimuthal anisotropy of the event results 
in large relative fluctuations of energy deposition in a ring 0. 

To conclude this section, let us discuss the influence of different factors on results for v 3 2 et 
determined with the present method. First of all, note that the theoretical absolute value of 
v 2 et and its 6-distribution are, of course, very mo del- dependent. For example, the effect of 
jet energy loss and corresponding azimuthal anisotropy are sensitive to the defined jet cone 
size [0, [11], |l^, [H|. Since we obtained here the energy loss only on the partonic level, for the 
finite jet cone size 8 ^ the effect should be less pronounced. Another reason for reducing jet 
azimuthal anisotropy can be sharp transverse collective expansion in events with large impact 
parameter values [] ||. However, the relative accuracy of v 2 et determination does not show any 
significant dependence on its absolute value. The reason for this is the following. The relative 
error for the method (|7]) using particle flow can be roughly estimated as a sum of two terms, 
which are proportional to (v J 2 e> 'N event )~ l and (v 2 art N^t)' 1 respectively (v 2 art is the coefficient 
of particle azimuthal anisotropy and N part is the mean multiplicity in the event). If the number 
of events is large enough, N event ^> 1, the main influence would be expected due to the second 
term, which does not depend on v 2 et . This is also the reason why the relative error becomes 
significant at low values of both the coefficient of particle azimuthal anisotropy and the particle 
multiplicity. The same conclusion will be valid for the method (|j) using energy flow, if the 
condition N part ^> N ring ^> 1 is fulfilled. 

Let us also note that in a real experimental situation, the pattern of azimuthal anisotropy 
can be more complicated due to non-flow correlations, finite accuracy of impact parameter 
determination, detector resolution effects, etc. On the other hand, one can try to improve ac- 
curacy of this technique considering, for example, the correlation between two equal multiplicity 
sub-events [13, or using results of work [|T9|. 



3 Note that the applicability of hydrodynamical model to very peripheral collisions is unclear. Moreover, the 
edge effects near the surface of the nucleus, impact parameter dependence of nuclear parton structure functions, 
early transverse expansion of the system and other potentially important phenomena for such collisions are 
beyond our consideration here. 

4 As it has been shown in recent work [|| , in the extreme scenario assuming instant transverse expansion of 
the medium, the geometrical anisotropy is strongly reduced, and the anisotropy of particles (at px > 2 — 10 
GeV/c under RHIC conditions) may drop below the observable level. In principle, this effect can have influence 
also on jet anisotropy at LHC conditions. Let us just mention that the possibility for realization of such 
scenario for transverse expansion is not obvious. Within original Bjorken approach [ p5| the rarefaction wave 
front moves with the sound velocity c s = \J dp/ de and the whole volume of the fluid appears to be involved in 
three-dimensional expansion in the time of order r ~ tq + Ra/c s . 
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4 Conclusions 



The strong interest is springing up to the azimuthal correlation measurements in ultrarelativistic 
heavy ion collisions. One of the main reasons is that the rescattering and energy loss of hard 
partons in azimuthally non-symmetric volume of dense quark-gluon matter can result in visible 
azimuthal anisotropy of high-p^ hadrons at RHIC and high-i?r jets at LHC. 

In jet case, the methodical advantage of azimuthal observables is that one needs to recon- 
struct only azimuthal position of jet without measuring total jet energy. One of the ways to 
perform the inclusive analysis of jet production as a function of azimuthal angle is event-by- 
event determination of the nuclear reaction plane angle. In the present paper we suggest the 
method for measurement of jet azimuthal anisotropy coefficients without reconstruction of the 
event plane. This technique is based on the calculation of correlations between the azimuthal 
position of jet axis and the angles of particles (not incorporated in the jet), azimuthal distri- 
bution of jets being described by the elliptic form. The method is generalized by introducing 
as weights the particle momenta or energy deposition in the calorimetric rings. In the latter 
case, we have illustrated the reliability of the present method in real physical situation under 
LHC conditions. The accuracy of the method improves with increasing multiplicity and parti- 
cle (energy) flow azimuthal anisotropy, and is practically independent of the absolute values of 
azimuthal anisotropy of the jet itself. 

To summarize, we believe that the present techniques may be useful for future data analysis 
in heavy ion collider experiments. 
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Figure 1: The impact parameter dependence of "theoretical" value of v 3 ^ including collisional 
and radiative energy loss (solid curve), and v J 2 et determined by the method (^) for dN ± /dy(y = 
0, b = 0) = 3000 (dotted curve) and 6000 (dashed curve). 
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